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Abstract 

According to the Quantum de Finetti Theorem, locally normal infi- 
nite particle states with Bose-Einstein symmetry can be represented as 
mixtures of infinite tensor powers of vector states. This note presents 
examples of infinite-particle states with Bose-Einstein symmetry that 
arise as limits of Gibbs ensembles on finite dimensional spaces, and 
displays their de Finetti representations. We consider Gibbs ensem- 
bles for systems of bosons in a finite dimensional setting and discover 
limits as the number of particles tends to infinity, provided the tem- 
perature is scaled in proportion to particle number. 



1 Introduction 

According to the Quantum de Finetti Theorem [2], locally normal infinite- 
particle states with Bose-Einstein symmetry can be represented as mixtures 
of infinite tensor powers of vector states. This note presents examples of 
infinite-particle states with Bose-Einstein symmetry that arise as limits of 
Gibbs ensembles on finite dimensional spaces, and displays their de Finetti 
representations. 

The central example is as follows. If the single-particle Hilbert space Ti 
is finite dimensional, the projector onto the symmetric subspace of the n- 
particle space can be normalized, and this defines the infinite-temperature 
ensemble for n bosons with single-particle space Ti. For each fixed m G N, the 
m-particle reduced density operators under the n-boson infinite-temperature 
ensembles converge, as n tends to infinity, to the density operator describing 
the m-particle statistics under a certain bosonic infinite-particle state ojq. 
The infinite-particle state ujq has a de Finetti representation as a mixture 
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of infinite tensor powers of vector states P Y , where v is a unit vector and 
P v = |v)(v| denotes the projector onto the span of v. In the de Finetti 
mixture for uq, the weight of the tensor power state p v ®°° is the probability 
density for v 6 Ti = C d+1 to equal 

v(p,0) = (e t9a VP~0i e^v^T,..., e**y/pi) (1) 

when p = (po,Pi, ■ ■ ■ ,Pd) is sampled uniformly from the (i-dimensional sim- 
plex A d and the phase angles 9i in 9 = (9 ,9 1 , . . . ,9a) are each sampled 
uniformly from [0,27r), independently of one another and of p. Thus the 
infinite-particle state uoq corresponds to the uniform probability measure on 
A, x [0,27r) d+1 . 

Similar limits are obtained for finite temperature Gibbs ensembles, pro- 
vided the temperature is scaled properly. Suppose H is a Hermitian operator 
on the single-particle space Ti = C d+1 and T n {(3) denotes the Gibbs canon- 
ical ensemble for n noninteracting bosons with single-particle Hamiltonian 
H at inverse temperature (3. Then, as n tends to infinity, the m-particle 
reduced density operators under T n ([3/n) converge to the m-particle den- 
sity of a certain bosonic infinite-particle state U p. The infinite-particle state 
up is an average of states P v ®°° with respect to the probability density on 
Ad x [0, 27r) d+1 that minimizes the "free energy" 

J J (v,Hv)f(p 1 9)dpd9 + - J j f(p,9)lnf(p,9)dpd9 , 

[0,270^+! d [0,2jr) d + 1 d 

where v = v(p, 9) is as in ([Tjl. We obtain similar results for bosons with "mean 
field" interactions, but again we must scale temperature in proportion to the 
number of particles. This stands in contrast to the analogous mean field 
limits for distinguishable particles, which are obtained without any peculiar 
scaling of temperature [Tj. 

The physical relevance of these facts is limited. On the one hand, they 
concern limits of canonical ensembles, which are appropriate when the num- 
ber of bosons is fixed, and therefore not appropriate for massless bosons (e.g., 
photons). On the other hand, massive bosons inhabit infinite dimensional 
Hilbert spaces, so to speak, whereas our results concern finite dimensional 
Hilbert spaces. However, the sort of ensemble we study is appropriate for 
(noninteracting) systems of n material bosons in thermal equilibrium, in case 
it is known that every one of these bosons is trapped in a potential well of 
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depth E. The statistical state of that system would be a conditional Gibbs 
ensemble, supported on the finite dimensional Hilbert space spanned by the 
symmetrized products of trapped (bound) states. Only noninteracting sys- 
tems of trapped bosons are considered, because the conditional Gibbs en- 
semble only makes sense if the Hamiltonian of the system commutes with 
the observable that every particle is trapped. 

Our results are presented in Section after a quick review of the Quan- 
tum de Finetti Theorem in the next section. 

2 The Quantum de Finetti Theorem 

Let TC be Hilbert space (which we will call the single-particle Hilbert space) 
and let TL® n denote the n-fold tensor power of TC (the n-particle Hilbert 
space). When 7r denotes a permutation of {1,2, ... ,n}, let U n denote the 
unitary "permutation" operator on Tt® 11 defined by 



For each n G N let D n be a density operator on the n-particle Hilbert space 
TC® n , the n-fold tensor power of TC. We want the density operators D n to be 
symmetric, and we assume 

(A) for all n, the density operator D n commutes with any permutation op- 
erator U n on H® n . 

We are especially interested here in systems of bosons, for which 

(B) for all n, D n U n = D n for any permutation operator U v on Ti® n . 

Condition (B) is stronger than (A). We also want the sequence {D n } of 
density operators to be consistent with respect to "subsampling" in the sense 
that 

(C) for all m < n, D n:m = D m , 

where D n:m denotes the m th partial trace of D n , i.e., the operator such that 



U w (x 1 <g> x 2 <E> ■ ■ • ® x n ) 



JCff(l) ® ^tt(2) ® • • • ® Xn(n) ■ 



n—m times 



Tr(D n:m A) 



Tr(D n (A®I (8)...® /)) 



for all A e B(Tt® m ). 
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The structure of sequences {D n } of density operators satisfying (C) and 
(A) or (B) is given by the quantum analogue of the de Finetti Theorem of 
probability theory [2|. Let p be a density operator on TC. A sequence {D n } 
of density operators of the form 

D 1 = p 
D 2 = p® p 

D 3 = p® p® p , et cetera (2) 

always satisfies (A) and (C), but it satisfies (B) and (C) if and only if p is a 
pure state, i.e., a rank one projector on TC. Roughly speaking, any sequence 
of density operators satisfying (A) and (C) is uniquely representable as a 
mixture of sequences of the form (j2J). That is, if {D n } satisfies (A) and (C) 
then there exists a unique probability measure p supported on the single- 
particle density operators such that 

D n = J p® n p(dp) (3) 

for all n. Furthermore, if {D n } satisfies (B) and (C), then the measure 
p(dp) in the integral representation © is even supported on the set of vector 
states p = P^p. This paraphrases the propositions of j2], ignoring the technical 
details; we now restate the results with more care. 

For m < n, let j mn denote the *-isomorphic embedding 

Jmn (B) = B®I® n ~ m 

of B(TC® m ) into B(TC® n ). The system of C* algebras B(H® n ) and isomorphic 
injections j mn has an inductive limit A. The inductive or direct limit in the 
category of C* algebras may be constructed as in |3J Proposition 11.4.1]. 
The inductive limit A is unique up to isomorphism, and for each n there is 
a ^isomorphism i n from B{7l® n ) into A such that i n jmn = i m for all m < n 
and the union of the images i n (B(H® n )) is dense in A. A sequence {D n } 
of density operators satisfying the conditions (C) can be used to define a 
continuous positive linear functional to on A by 

u(i n (B)) = Tr(D n B) MB e B(H® n ) . (4) 

This is well-defined thanks to the consistency conditions (C) and the density 
of Ui n (B(H® n )) in A. In particular, u(e) = 1, where e is the identity element 
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of the C* algebra A. If {D n } satisfies (A) as well as (C) then uj is symmetric 
in the sense that 

uiiniUirBU*)) = u{i n {B)) (5) 

for all n, all B G B(H,® n ), and all it G IT n , the set of permutations of 
{1, 2, . . . , n}. The set of all "symmetric states" on A, i.e., the set 

SS = | a; G A* u(e) = 1 and u(x*x) > Va; G A and uj satisfies (0) j , 

is a convex subset of the Banach dual A* of A, and it is compact with respect 
to the weak* topology. Let SSi denote the space of single-particle states, 
i.e., the set 

SS! = jp G B{H)* p(I) = 1 and uj(A*A) > VA G B(H)^ 

endowed with the relative weak* topology it inherits as a subset of the Banach 
dual B(7i)* of B(TC) . It was first shown in jl] that each uj G SS has a unique 
representation as an integral of product states 



uj = / p ® p ® p ® fi(dp) = / p®°°p(rfp) , (6) 

JsSi JSS! 

where p is a probability measure on the a-algebra T\ generated by the in- 
tersections with SS\ of weak* open sets in B(TL)* . We sketch a proof of this, 
following reference |2]: First, the extreme points of SS are identified as the 
product states p®°° . Thus, the set of extreme points is the image of the com- 
pact space SSi under the continuous injection p i — > p®°°, and it follows that 
the extreme set is closed in SS. The existence of an integral representation 
(JBJ) is then a consequence of the Krein-Milman Theorem, and its uniqueness 
is shown in [2] by a direct argument. 

It is further shown in [2] that the measure p(dp) appearing in the integral 
representation © of uj is supported on the measurable subset of normal 
states on B{7i) if uj is determined, as in formula (J1J above, by sequences 
of density operators satisfying (A) and (C). If, in addition, the sequence of 
density operators defining uj satisfies (B), then the measure p(dp) is even 
supported on the vector states p(A) = (ip,Aip) with ||?/>|| = 1. 



3 Examples of bosonic de Finetti states 

In this section we exhibit some sequences {D n } satisfying (B) and (C) that 
are obtained from natural statistical ensembles. In all of these examples, the 
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single-particle Hilbert space Ti. is finite dimensional. After introducing the 
notation, we will state all of our results before proceeding to their proofs. 

Let Ti = C d+1 and let Ti^ n > denote the subspace of symmetric vectors in 
Ti® 11 . Let S n denote the symmetrizing projector 

= 1 £ (7) 

7ren„ 

from Ti® n onto Ti {n) . We now introduce notation for the occupation num- 
ber basis of TiS n ' relative to a fixed orthonormal (ordered) basis {ej} of Ti. 
Let n = (no, Tii, . . . , no) be an ordered d + 1-tuple of nonnegative integers 
(occupation numbers) and let #n denote Yl n j- We use the notation 

I) = «vn«<! 

/ i=0 

for multinomial coefficients. The vector 
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* n = J ( U ) S n (er° ® ef ni ® • • - ® ef*) ( 8 ) 



is a unit vector in Ti^ n \ and the set of vectors {^n | #n = n} is an orthonor- 
mal basis of H.( n \ Let P n denote the rank-one projector onto the span of \l/ n : 

P n $ = (9) 

We begin by considering the "uniformly mixed" density operators sup- 
ported on Tiy^\ 

Proposition 1 Let S n denote the symmetrizing projector |7L). For each m, 

s —i^Tvk E -= e {(™)/ A n™ P )}p m , do) 

n m:#m=m v 7 Jfll i=0 

where Xd(dp) denotes normalized Lebesgue measure on the d-dimensional sim- 
plex 



A 



a 

d = = (Po,Pi, • • • ,Pd) e JR d+1 | < pi i = l,2,...,d and ^Pi = 1 

i=0 
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The sequence {S m } satisfies (B) and (C) of Sectional By the Quantum de 
Finetti Theorem, there exists a measure /i supported on the pure states on 
C d+1 such that 

S m = I P® m fi(dP) 



for all m G N. This measure can be described as follows. Define the map 

v : A d x [0, 2n) d+1 — ► C d+1 

by 

d 

v(po,Pi,...,Pd,0 ,6 1 ,...,6 d ) = ^e^y/Pj ej (H) 

3=0 

where {e{\ is the standard basis of C d+l . The map v is many-one onto the set 
of unit vectors in C d+1 . The probability measure p,(dP) is the one induced 
via v from the uniform measure 

on Ad x [0, 27r) d+1 . In other words, 
Proposition 2 The density operator hlf^l equals 
c times 

/ {P V{P ,8) P vM )a{dO)X d {dp) . (12) 

A d J[0,2Tr) d + 1 

Next we consider Gibbs ensembles for noninter acting systems of bosons. 

Let 

n 

H n = ^T, (13) 
1=1 

be the Hamiltonian for n noninteracting bosons with single-particle space 
Ji = C d+1 . Let {ej} be an orthonormal basis of 7i consisting of eigenvectors 
of the single-particle operator T, so that Tej = e^j. The Gibbs density 
operator for the n boson system is 

d d 

r„(/3) = — H e ^ nmp * with = E U e ~ Pn%tl ■ 

n ^ n:#n=n i=0 n:#n=n i=0 

(14) 

An interesting limit is attained if temperature is scaled in proportion to n 
as n — > oo. If the temperature is not scaled as n — > oo then a sort of 
Bose-Einstein condensation is attained in the limit. 
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Proposition 3 Let H n be the noninteracting Hamiltonian hlty and let T n ((3) 
denote the Gibbs density Let {ej} be an orthonormal basis ofH con- 

sisting of eigenvectors of the single-particle operator T, so that Tej = ejCj. 

(i) For each m 6 N, the limit lim T n . m (f3/n) exists and equals 

n ->oc 

£ {(") / flpT'Z^fle-^X d (dp)}p m 

m:#m=m ^ ' J A d i=0 1=0 

with Zp 1 = J Ad nto ex P(-^ e ^)^d( rf P)- 

(ii) If e < e 1 < ■ ■ ■ < e d , then for each m GN, 



lim T n ., m (f3) = P( m , ,...,o) = P ( 



e 

Finally, we consider systems with two-particle interactions in the "mean 
field" scaling. Let V be a Hamiltonian operator on 7i ® 7i such that V(x <8> 
y) = ^(y ® x ) f° r an x,y ETC. For n > 2, define the Hamiltonian 

n 

H n = J> + — - £ V«, (15) 

where V^- denotes the operator obtained by applying V to the 2 th and j th 
factors of TL® n . For any n G N and any ,3 6 1, the n-particle Gibbs density 
at inverse temperature j3 for the Hamiltonian (|15|) is 

Proposition 4 Let r n (/3) denote the Gibbs density Mb}) . For each m, the 
limit 

G m = lim {r n (/5/n)}. m 

exists and defines a density operator on (C d+1 ) lg>n . The de Finetti represen- 
tation of G m is 



-L / / P v ® ■ ■ ■ ® P v e^ TP ^ v ^ P ^a(d6)X d (dp) 

Z P JA d J[0,2ir) d + 1 

with v = v(p, 0) as m iii)) and 

Z p = f [ e-^ TP ^ v ^ p ^a(d9)\ d (dp) . 



8 



The rest of this section is devoted to proving the above propositions. 
Recall the definition (jUJ) of the projectors P n . For each n & N, let p n be 
an n-particle density 

Pn = 22 W n( n )Pn , 
n:#n=n 

where w n is a probability measure on the set {n|#n = n}. Each probability 
measure w n can be associated with the discrete probability measure 



L) r , 



n:#n=ra 



on the d-dimensional simplex A^. It may be verified that 

-1 d 



P 

n -r 



e n 

m:#m=m i=0 



m. 



(this equals if any m 8 > rii for any i), and therefore 



E I. E n 

m:#m=m n:#n=n 



8=0 



P, 



E 

m:#m=m 



rrd th/th _ 1 

w ^ n > i(i-l)(i- 

n:#n=n «- yv 



> n n 



m— 1 ' 



(17) 



The coefficient of P m in (|T7jl may be written 

(m) Sa f n ^ Wn ( dp } ' 

where 



mi— 1 ' 



/n(p) — H{pi>(mj-l)/n Vi}(p) jw^ 2 \ q m-l \ ■ 

The functions / n (p) converge uniformly to n^Lo^r 4 011 A<j. Therefore, if 
converges weakly to some probability measure u(dp) on Ad, then 



Mr, 



lim p n: . 



E 

m:#:m=m 



Y[pTMdp) Pm • 



;is) 



i=0 



The probability measures on corresponding to the Gibbs density op- 
erators (|14|) for noninteracting bosons are 

d 

n:#n=n i=0 

If all of the eigenvalues of T are equal, then the measures ()19|) converge weakly 
to Xd(dp), the uniform probability measure on the simplex, but if eo is strictly 
smaller than all of the other eigenvalues of T, then the measures (j!9|) converge 
weakly to S(p — (1, 0, . . . , 0)), a point-mass at the lowest energy vertex of 
the simplex. This convergence implies Propositions Q and assertion (ii) of 
Proposition El by formula ([18)1 . On the other hand, the probability measures 
corresponding to the Gibbs density operators T n (f3/n) for noninteracting 
bosons are 

d 

w « = z i E U e ~^ /n6 (p-^ > 

n:#:n=n i=0 

and these converge weakly to 

d 

with Z/3 = f A Yii=o ex v(— P^iPi) ^d(dp) . This proves assertion (i) of Proposi- 
tion El 

To prove Proposition El we will show that ()12|) and (|10|) are equal. Define 
the rank-one operators Qjk{x) = {&ki x ) e j- From ([lip. 

d 

j,k=0 

and therefore Pv{p,6)® m equals 

d d d 

E EII VPj^ ei{9jr ~ dhr) Qhk, ® Q J2 fc 2 ® • • • ® %^ . (20) 

ill— Jm=0 fci,...,fe m =0 r=0 

For i = 0, 1, . . . , d, let N { : {0, 1, . . . , d} m — ► N be defined by 

Ni(x u x 2 , ■ ■ ■ ,x m ) = #{k e {1,2, ...,m} : x fc =i} 
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and define 

N(x 1 ,x 2 , . . . ,x m ) = (N Q (xi,x 2 , ■ ■ ■ ,x m ), ■ ■ ■ ,N d (xi,x 2 , . . . ,x m )) 
If N(ji, ...,j m ) = N(ki, ...,k m ) then 



[0,2^+1 r=Q 



JJe^-^aide) = 1 , 

r=0 

but otherwise it equals 0. Thus, from (|20|). 

P v(p ,ef m a(d9)\ d (dp) 

P m 



m:#m=m d i=0 

N(ji,...,j m )=m N(k 1 ,...,k m )=m 



= e / 4 n^*)(™V. 

m:#m=m J A <* i=0 7 



by the definition ® of P m . This proves Proposition El 

Finally, we derive Proposition 0] from Proposition El Define W = T®/ + 
/ (g) T + V. Then the Hamiltonian ()15|) can be written 



zrr E «V 



#n = 

n 

l<i<j<n 

We claim that 

lim — — {{H n yZ n } :m = 2~ 3 {W m+ i >m+2 W m+3 ,m+4 ' ' ' W m j r2 j-i m j r2 jS m j r2 j\ 

(21) 

for each j, m G N. This is so because (H n Y contains (j£) J terms of the form 
(n — l) _: 'W 01 6 1 K / a2 6 2 • • - Wo,;,-, and, when n is large, the majority of these 
terms are such that the indices a 1; & 1( . . . , a,,-, are all distinct and greater 
than m. The sum of the remaining terms in (H n Y is o(n J ) and does not 
contribute to the limit (|21j). By the symmetry of E n , 

{W ai biW a2 b 2 ■ ■ ■ Wajbj^n} . m = {W / m+l,m+2 ■ • ■ W / m +2j-l,m+2jE n } :m 

= {W / m+l,m+2 ' ' ' W / m+2j-l,m+2jEn:m+2j }. m 
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if ai, bi, . . . , cij, bj are all distinct and greater than m. There are asymptoti- 
cally n 2j /2 such terms, so (|21j) follows from Proposition ^ 
Now, to prove Proposition 0J expand 



(22) 



1 00 1 1 

and take the m" 1 partial trace: 

^{^-'^}„ = gjf(-«'»- , {(^^} m - 

The term of the series in ()22j) converges to 

( — l) J ~y^j {^+1^+2^+3,01+4 ' ' 'W rn+ 2j-l : m+2jS rn +2j}. m 

as n — ► oo by (}2*T]) and is bounded by 4y /3 J \\W |p uniformly in n. Since 
the series in f!22|) are majorized by the convergent series Ylj II ^IP an d 
converge term-by-term as n — > oo, it follows that 



n i=o J ' 

(23) 

Substituting the integral representations (JT2J) for S m+ 2j into (|2lj|l yields 



l im 1 {e-^^-sJ 
n-*oo TrL n 



j=0 ^ 17 A d J iO^) d+1 



,1. 



e -^2 r& ( WP v( P , e )®^( P , e ))p v(p)e) ®m .(^) A ^ fip ) _ 

A d J[0,2n) d + 1 

Proposition 0] follows from the preceding equation and the definition (jlfij) of 

r n (/3). 
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